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Abstract. The Manley-Rowe constants of motion (MRC) are conservation laws writ- 
ten out for a dynamical system describing the time evolution of the amplitudes in res- 
onant triad. In this paper we extend the concept of MRC to resonance clusters of any 
form yielding generalized Manley-Rowe constants (gMRC) and give a constructive 
method how to compute them. We also give details of a Mathematica implementation 
of this method. While MRC provide integrability of the underlying dynamical sys- 
tem, gMRC generally do not but may be used for qualitative and numerical study of 
dynamical systems describing generic resonance clusters. 
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1 Introduction 

Weakly nonlinear wave systems showing resonances can be found in a great variety 
of physical systems, among them surface water waves, atmospheric planetary waves, 
plasma drift waves, etc. In Hamiltonian formulation the equation of motion in Fourier 
space can be written out as 

iB k = dH/dBl (1.1) 

where is the amplitude of the Fourier mode corresponding to wavevector k and the 
Hamiltonian % is represented as an expansion in powers of terms Hp each representing 
all products of amplitudes B^: 

H = H2+H int , Hi nt =n 3 +H 4 +n 5 +.... (1.2) 

The quadratic Hamiltonian H2 describes linear motion and the interaction Hamiltonian 
Hint describes nonlinear interaction of waves. 
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If the cubic Hamiltonian 7i 3 ^0, Hm^Hs, three-wave interaction is dominant and the 
main contribution to the nonlinear evolution comes from triads of waves each satisfying 
the resonance conditions 

w(k 1 )+&;(k2)-a;(k3)=0 / k!+k 2 -k 3 = 0. (1.3) 

The corresponding dynamical system describing time evolution of the amplitudes Bj of 
a triad reads 

Bj = ZB* 2 B 3 , B 2 = ZB\B 3 , B 3 = -ZB X B 2 , (1.4) 

the interaction coefficient Z 7^ being determined from a known function over the solu- 
tions of dl.3b . 

System (|1.4[) possesses three conservation laws, first found in [1J called Manley-Rowe 
constants (MRC) of motion. For a resonant triad, the MRC may be written as: 

fl3= |£>i| 2 +|B3| 2 / ^23 = I S2 1 2 -h |f>3| 2 / hi = |Bl| 2— |£>2| 2 / (1.5) 

They are linearly dependent, but any two of them form a linearly independent subset, 
enough to provide the integrability of (|1.4|) , |2|. 

In a 3-wave system, a resonant triad is called a primary resonance cluster, and clusters 
consisting of more than 3 modes (common clusters) may be decomposed into triads hav- 
ing joint modes, [3]. Common clusters describe many real physical phenomena: in non- 
linear water wave systems |4J, in piezoelectric semiconductors |5j, in geophysics |6tlT0 | , 
in optics |7| etc. 

In this paper we present a simple constructive method for deducing generalized MRC 
(gMRC) for a common resonance cluster. The gMRC are polynomials on \Bj\ 2 with in- 
teger coefficients. The method has been developed by EK for the course on nonlinear 
resonance analysis held at the J. Kepler University, Linz, since 2005. A Mathematica code 
based on this method has been written by Loredana Tec and can be downloaded from [8 J. 
The method has been used afterwards by the author, collaborators and students but its 
description has not been published. Recent publication [18J where gMRC are studied by 
far more complicated methods indicates that a simpler algorithm and available program 
code may be of interest for physicists working in the area of discrete wave turbulence. 

2 Computation of gMRC 

2.1 Computation of MRC 

Let us rewrite (|1.4|) as 



B\ = ZB\ B 3 , B 2 = ZB\ B3, B 3 = —ZB\ B 2 , 

B- * 7B D* jS* 7D D* B* 713*13* 
j — Z J i52"3, t->2 — i-iD\D 3l l> 3 — — ZjD^ 



(2.1) 
(2.2) 
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where the second equation is the complex conjugate of the first. As BjBJ = \ Bj\ 2 and 

d(B ; B* )/dt = BjdB* /dt+BjdBj/dt, (2.3) 
multiplication of each equation for Bj and BJ by Bj and By, respectively yields 

{d|B 1 | 2 /df = ZB?B|B 3 + ZB 1 B 2 B* 
d|B 2 | 2 /df = ZB*B 2 *B 3 + ZB 1 B 2 B^, 
d|B 3 | 2 /di=-ZBiB 2 B^-ZB*B|B 3 . 



'd|Bi| 2 /df-d|B 2 | 2 /dt = 0, 
d|Bi| 2 /dt + d|B 3 | 2 /dt = 0, 
i d|B 2 | 2 /df + d|B 3 | 2 /dt = 0, 



7 12 = |B 1 | 2 -|B 2 | 2 : 

J 13 =|B 1 | 2 +|B 3 | 2 : 
U 23 =|B 2 | 2 +|B 3 | 2 : 



: const, 
= const, 
- const, 



(2.4) 



which is the set of MRC found originally in [1J. 



2.2 Computation of gMRC for a two triad cluster 

Regard a simple common cluster in a three-wave system formed by two triads a and b 
connected by the low frequency modes B\ a = B\i. The corresponding dynamical system 
reads 



B\ a — ZaB'^Bsa + ZijB^B^i,, 

B'* 7 U D* i 7 u P* 
1 fl — Z a #2fl + Z 6 2fo # 3ft / 



U 7 D* U 72* 7 U 12* 

7 T> U D* 7 U* D* 

3a — —^aD\aD2a, £> 3fl — — Z fl B lfl £> 2fl , 



u 7 D* K P* — 

D2b — ^b0 la D3b, D 2b - 
B3b = —Zj,Bi a B 2 i,,Bj h 



ZfoBJ fl B 3 b, 

ZD* D* 



Considerations similar to those above yield three gMRC: 



(2.5) 



^d|B 3b | 2 /d£: 

d(|B lfl | 2 + |B 3fl | 2 + |B 3fc | 2 )/d£ = 0, 
d(|B 2fl | 2 +|B 3fl | 2 )/d£ = 0, 



d|Bia| 2 / dt = Z a B\ a B\ a B3 a + Z} ) B\ a B1 lh By, 
+Z a B\ a B 2a Bg H + Z(, B\ a B 2 b B| b , 
d | B 2h 1 2 / d£ = Z a B\ a B 2a B 3fl + Z a B\ a B 2a B^, 
d|B 3fl | 2 / dt = —Z a Bi a B 2a B^ a — Z a B\ a B\ a B3 a . 
d | B 2b | 2 / dt = Z b B\ a B| b B 3b + Z h Bi a B 2b B^ b , 

ZD* D* D 



|d(|B 



ZbBi a B 2 bB^ b 

f|Bi„| 2 + |B 3fl | 2 
|B 2fl | 2 + |B 3fl | 2 



|B 3b | 2 )/d£ = 



3 2b| 



IB 



3b 



|B 3b | 2 = const, 
const, 
2 = const. 



(2.6) 
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It is important to note that the choice of the connecting mode influences the form of the 
generalized MRCs. Indeed, for the connection B\ a = B^, we got the set of constants (|2.6[) 
whereas the connection B\ a = B^ yields a different set of MRCs: 

( J 

< hi\b 

k I\a,b 

The sets (12.6ft and (I2.7D are linearly independent and the corresponding dynamical sys- 
tems are not equivalent. 

As it was first shown in [9 1 the connections within a cluster of triads may be classified 
into three connection types depending on whether a connection involves both or one or 
none of the high frequency modes B 3a , B 3 /,. A high frequency mode is unstable due 
to the criterion of decay instability and is called active or A-mode while low frequency 
modes £>i fl ,£>2a/£>ifc/£>2fa are called passive or P-modes. Each connection within a common 
cluster can be either of AA-, AP- or PP-type (see |2| for more details). In the NR-diagram 
representation they are shown as two bold, bold-dashed and dashed-dashed half edges 
correspondingly and each triad is shown as a triangle 

Examples of all possible NR-diagrams for 2-triad clusters are shown in FigHJ For the 
NR-diagrams shown in the upper panel on the left and in the middle the sets of MRCs are 
given in (|2.6|) and \2.7\ respectively. The variety of 3-triad clusters is substantially richer 
(a few dozen). Examples of NR-diagrams for some possible clusters are shown in Figj2j 
In real physical applications, clusters of 5 to 10 triads and more are common, e.g. [9-12J. 
The sheer voluminous of computations for deducing MRCs makes the use of a computer 
algebra system mandatory. 

Another important issue is that as any linear combination of gMRCs is a gMRC, gM- 
RCs are not unique; any of the constants in (|2.6|) can be replaced by some of their linear 
combinations: 

; 1 = |Bi (I | 2 -|B 2 «| 2 + |B 3 b| 2 ,or (2.8) 
h = I Bia 1 2 + 1 B 2a | 2 + 2 1 B 3a | 2 + 1 B 3b | 2 , or else. (2.9) 

However, for qualitative and numerical studies we are naturally interested in computing 
"a minimal set" of gMRC, i.e a set of the shortest possible form which contains specific 



l B 2|«| 2 + l B 3|fl| 2 ' 

l B l|b| 2_ l B 2|b| 2 / ( 2 - 7 ) 
B 3\a\ +l B 3|fcl • 
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Figure 2: NR-diagrams for some 3-triad clusters 



variables. For instance, for the dynamical system (|2.5|) , the set 

I h3\a = \^2a\ 2 +\B3a\ 2 , ^316 = I ^2b | 2 + I ^3b | 2 / 



I\ab— |£>lfl| 2 +|f>3a 



IB 



3^ 



(2.10) 



is minimal while the set 



{h3\a— \Bla\ 2 + \B3a\ 2 , Il3\b ~ I ^2b \ 2 + I ^3b | 2 / 
h = |Bifl| 2 +|B2fl| 2 +2|B3 fl | 2 +|B 3 i,| 2 . 

is not. Minimal sets of gMRC can be computed using Grobner bases. 



(2.11) 



3 Implementation details 

The computation of gMRC has been implemented in Mathematica as a package NonlinearResonance 
including the following procedures (algorithms are given below): 

• DynamicalSystem for computing dynamical systems for common clusters in a three 
wave system; 

• ConservationLawsLinComb for computing the MRCs for a given cluster by consid- 
ering all linear combinations of the variables from the corresponding dynamical 
system having coefficients from {0,±1}. The set of conservation laws is obtained 
by computing the null space of the set of these vertices; 



• ConservationLawsGb for computing a Grobner basis of the ideal of the MRCs for a 
common cluster. 
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Algorithm 1 DynamicalSystem 

Require: L = ((L//)y_ B );_i a list of triads representing a cluster. 
Ensure: The dynamical system R of the cluster. 

1: M:= (ZiB^ b ,Bi c ,ZiB^ a Bi c ,—ZiBi a B lb ,...,Z s B* b ,B sc ,Z s B* a B sc ,—Z s B sa B sb ); 
N:= (Bi a ,Bu,,B-[ c ,...,B S a,B s i,,B S c); 
L'-= (Li a ,Lu,,Li c ,...,L sa ,L s ),,L S c); 
2: for i from 3s by — 1 to 2 do 
3: for from i — 1 by — 1 to 1 do 
4: if Li = Lj then 

5: Mj-.= Mi+Mf, 

Mi-.= NULL; 

In M substitute N; by N y ; 

Ni-. = NULL; 

;:=l; 

6: end if 

7: end for 
8: end for 

9: for i from 1 to |M| do 
10: Ri-.^N^Mi); 

R i+m ~(f t N*=M*y, 
11: end for 
12: return R; 



Algorithm 2 ConservationLawsLinComb 



Require: L a list of triads representing a cluster. 
Ensure: A minimal set M of MRCs for the cluster. 
1: L:=DynamicalSystem(L); 

2: Let n be the number of nodes appearing in the dynamical system L, say B\ lt ...,B\ ; 
3: forl<i<wdo 

4: |(I B A,| 2 ):=Ba,|(B:,) + |(B A ,)B:,, where |(B A ,),|(B* i ) are determined by the dynamical system; 
5: end for 

6: Determine a basis M of 

{(^...^MCiir | j t (v 1 \B A f + ...+v n \B A f)=0y, 



7: return M; 
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Algorithm 3 ConservationLawsGb 
Require: L a list of triads representing a cluster. 

Ensure: A Grobner basis G of the ideal of conservation laws and the corresponding term order 
1: L:=DynamicalSystem(L); 

2: Let n be the number of nodes appearing in the dynamical system L, say Bx 1 ,.-./ B X n 'i 
3: f or 1 < i < n do 

4: id B A,| 2 ) ~ B Xift( B V + Tt( B k) B li> where Tt( B *,)'Tt( B V are determined by the dynamical system; 
5: end for 

6: Choose a term order -< such that for t x 6 [f t (\ B M | 2 ),- - ^ ( | Ba„ | 2 )], 1 7^ *2 6 [B A] ,... / B An ,B* i ,... / B*J, t 3 e 
[^(|B Al | 2 ),... / ^(|B A J 2 ),B Al ,... / B Au ,B* i ,...,B A J we have £ a -< t 2 *3 and such that -< restricted to 
li(\ B M\ 2 ),-4(\ B X n \ 2 )}is degree lexicographic; 
7: Compute a Grobner basis G of the ideal <^(|B Ai | 2 )-B A .^(B^) + ^(B A .)B* i | i = l,...,n> < 

Q[|(|BaJ 2 ) |(|B A J 2 ),B A]/ ... / B An/ B Ai/ ...,B*Jwithrespectto^; 

8: G = GnQ[|(|B Al | 2 ),... / |(|B A „| 2 )]; 
9: Remove all elements of G with degree > 2; 
10: In G substitute ^(|B A; | 2 ) by |B Ai | 2 for i = l,...,n; 
11: return Gand^;; 



Each of those takes as input a list of triads, where every triad is represented by a list of 
its vertices. The procedure DynamicalSystem returns the dynamical system of the cluster 
specified by the input. The procedure ConservationLawsLinComb returns a "minimal" 
basis of the set of conservation laws - i.e., a set of MRCs - corresponding to the input 
cluster. The procedure ConservationLawsGb returns a normalized reduced Grobner basis 
of the ideal of conservation laws with respect to the lexicographic term ordering on the 
\Bjj\ 2 s such that 

l^| 2 >l%| 2 { ifZ< "° r 

^if/ = fcand;'</ 

This is usually computed faster than the basis returned by conservationlaws but its ele- 
ments need not only have coefficients ±1 and 0. Computation time for arbitrary clusters 
consisting of 20-30 triads is about a few seconds. 

4 Discussion 

The construction of any constants of motion aims usually at establishing the integrability 
of the underlying dynamical system. For instance, integrability of a two triad cluster 
with £>3 n = £>3(, and Z a /Zj, = 2 is proven in [13] where an additional constant of motion 
(a polynomial of sixth order) has been constructed via direct search based on irreducible 
forms. On the other hand, a two triad cluster with B\ a = B\b and Z fl /Zfc = 3/4 shows 
chaotic energy exchange among the modes of the cluster (see Poincare sections for this 
case in |2J, Fig.4.6). 

Examples of integrable dynamical systems are known only for resonance clusters of 
a very special form and/or for special initial conditions. Moreover, there is no general 
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method to establish integrability of a dynamical system and usually only a small set of 
constants of motion of any form (other than gMRC) can be found, e.g. 1121 (13141711 . 

The algorithm described in this paper allows to produce in a matter of seconds (mak- 
ing use of a computer algebra system) a set of gMRC which is appropriate for qualitative 
study of the energy transport in a resonance cluster, for simplifying numerical study by 
reducing the number of variables in a dynamical system or for checking stability of the 
chosen numerical scheme. 

Generalized MRCs can be analogously constructed for a 4-wave system where the 
primary clusters are quartets. The explicit form of the gMRCs for a quartet can be found 
in 0. 
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